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Fermionic Loops in Numerical Stochastic Perturbation Theory 
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We discuss the inclusion of fermionic loops contributions in Numerical Stochastic Perturbation Theory for 
Lattice Gauge Theories. We show how the algorithm implementation is in principle straightforward and report 
on the status of the project. 



1. Introduction 

In recent years the Numerical implementation 
of Stochastic Perturbation Theory (NSPT) was 
introduced, which was able to reach unprece- 
dented high orders in perturbative expansions in 
Lattice Gauge Theories (LGT) . Till now the main 
limitation of the method has been the quenched 
approximation, an inclusion of fermionic loops 
contributions missing. We can now fill this gap. 
Let us first of all remind the basics of the method. 
NSPT § comes (almost for free . . . ) as an appli- 
cation of the Stochastic Quantization scheme of 
Parisi and Wu [|] . We briefly sketch it in the con- 
text of LGT. The Cornell group showed ||| that 
a Langevin equation for LGT can be formulated 
as 



U^x-T + e)=e- F V^^U tl {x-T) 



(1) 



where 
F[U(T),r,) 



i 

^T\eV^S[U] + V~^)- (2) 



In the previous formulae S[U] is the action for 
the gauge fields, e is a time step and t] is a gaus- 
sian noise, while V^. „ is a Lie derivative on the 
group, whose definition can be easily understood 
in terms of 

f(e a - T U) = f(U) + a l V l f{U) + 0(a 2 ). (3) 

Eq. (Q) should be understood as follows. One has 
at hand the theory described by the action S[U] 



and the goal is to compute observables in terms 
of the functional integral, in which the measure 
is dictated by exp(— S[U]) (we adhere to an eu- 
clidean formulation). In Eq. ([!]) an extra depen- 
dence on a new parameter r is imposed on the 
fields. This parameter can be thought of as a 
stochastic time in which an evolution takes places 
according to Eq. ([!]), which is stochastic due to 
the presence of rj. Now the key points are the 
gaussian nature of r\ and the fact that the so called 
drift term in the F appearing in Eq. ([!]) is given 
by the equations of motion. One can show that 
because of that the Langevin equation describes 
a stochastic process whose (asymptotic) equilib- 
rium distribution is given just by the measure one 
is interested in (exp(— S[U])). Therefore one can 
trade expectation values with respect to the lat- 
ter for means over the stochastic time evolution 
dictated by Eq. (0). In the case of pure gauge 
Wilson action 



2N 



(4) 



Eq. (|) reads 



AN 

i UpBU^(x) 

(the subscript tr asks for the traceless part). One 
can recognize in Eq. (||) a pretty local nature, 
as expected. It is well known [|| that one can 
get a Stochastic Perturbation Theory from the 
Langevin equation. NSPT gets it by directly ex- 
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panding the fields as 



U^(x) = expL4 M (a 



exp 



E^~ fc/2 4 fc) (-) 



.fe>0 



(6) 



l + ^/3- fc / 2 C/W(x). 



fe>0 

Such an expansion has to be thought of as a for- 
mal series motivated by the dependence of the 
solution to Eq. (|l|) on the coupling {(3~ k / 2 ~ g). 
As a consequence every quantity depending on 
the fields can now be given a similar power ex- 
pansion such as (for example) 



f[u] = j2r k/2 F {k) . 



(7) 



k>Q 



Eq. (|l|) gets translated into a hierarchy of equa- 
tions exactly truncable at any given order and 
suitable for a numerical integration on a com- 
puter. Each order of Perturbation Theory is now 
obtained from means over stochastic time evo- 
lution of suitable (maybe complicate) composite 
operators (just like those appearing in Eq. (fj])). 
Notice that having decompactified the formula- 
tion in Eq. (^), divergencies in the non-gauge- 
invariant sector show up, which are cured by 
Stochastic Gauge Fixing. 

2. Introducing fermions 

In principle introducing fermions is straightfor- 
ward. One needs to face a new measure in the 
functional integral of the form 



e- s ° d C tM = ( r^- TrlnM \ 



(8) 



We write M for the fermion matrix to recall that 
other actions have the same structure, for exam- 
ple the Faddeev — Popov action Q . One would in 
principle simply need to replace 



V^S G -V^Tr In AT = 

V* S G -Tr((V* M)M- X ). (9) 



As a matter of fact, one has to face an inverse, i.e. 
non-locality. The Cornell group S gave to the 



problem a solution which was in a sense a proto- 
type for fcrmionic simulations. This amounts to 
rewrite in Eq. (H) 



(10) 



in which a summation over repeated indices is to 
be understood. One should also keep in mind that 
k,l,n arc mufti indices. A new gaussian field has 
been introduced, normalized as (£<£j) = Sij. The 
evolution of the process will now also average over 
£, resulting in 



= [K,„S G - Tr ((V^M)M- 1 )] 
= V* „ [So ~ Tr (InM)] 



(11) 



which is exactly what one is interested in. Life 
is now pretty easier, since one simply needs to 
invert the fermion matrix on a source solving the 
system 

M kl ibi = & (12) 
in terms of whose solution the evolution is local: 



Fi = 



Re&^V* M)mM1 • (I 3 ) 



It takes a little time to realize that life is in a 
sense even easier in NSPT. To understand why, 
one should first of all remember that also the 
fcrmionic matrix gets expanded as a power series 



fe>0 



(14) 



The inverse of such a matrix is easy to compute 
M" 1 = ^/3- fe / 2 M- lW 



k=0 



J2$- k/2 M- m . (15) 

fc>0 



The notation enlightens the fact "the zeroth- 
order of the inverse is the inverse of the zeroth- 
order" . Other orders are not difficult to compute, 
resulting out of a simple recursion 



M 
M 



= -M^M^M^ 1 
-i (2) = -M^M^M^ 1 



(16) 
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-m(°) _1 mWm- i(1) 
-mW^mWm- 1 ™ 



An NSPT algorithm for fermions is now easy to 
implement in force of the four considerations we 
now proceed to make. 

2.1. Locality 

First of all, from the new random fields £ define 
the (power-expanded) fields 

^=A/- l(j) £ (17) 

in terms of which the (I + j)th contribution to 
Eq. (0) reads 

£kW, M M)®tP+{l~j, ifl^j) (18) 

The expression of (y i x fi M)fj is easy to work out. 
In order not to obscure the main point with trivial 
algebra, we do not write it down, only stressing 
that it is as local as its non-perturbative coun- 
terpart, as obviously expected. 

2.2. An easy recursive formula 

Having made the point that a simple recursive 
formula holds for Eq. (jl^) , it is straightforward to 
notice that a similar relation holds for the com- 
putation of the fields 



£(0) 








(19) 


fd) 








£(2) 




M (2)£(0) _ 






£(3) 











The message from Eq. (|l^) is quite clear: at ev- 
ery order only one proper inverse is needed and 
everything that is left comes out of a simple recur- 
sive relation in terms of already computed (local) 
quantities. On top of that, the main point still 
needs to be made as far as the matrix is 
concerned. 



2.3. The call for an FFT 

Being the Oth order, M^ ' of course does 
not depend on the fields: it is actually the stan- 
dard Feynmann propagator, which is diagonal in 
Fourier space. Notice that the latter observa- 
tion holds for any fermionic matrix and also for 
the case of the Faddeev-Popov matrix With 
this respect the method is quite general, opening 
the way to multiple applications, among which 
Neuberger fermions, whose perturbative expan- 
sions have till now only been pioneered. Given 
the above observations, the obvious way to im- 
plement the construction of the fermionic contri- 
bution to the drift F is to go back and forth from 
Fourier space. This of course calls for an effi- 
cient FFT. The implementation we are working 
on is the 4-d version of the algorithm described 
in || . In general, one can always think of a multi- 
dimensional FFT as the result of subsequent ap- 
plications of 1-d FFT. This approach is the ob- 
vious choice on an architecture such that of the 
APE family computers, in which case the algo- 
rithm basically amounts to a 1-d (local) FFT in- 
terlaced with transpositions^]. Implementing an 
FFT in the NSPT context also opens the way to 
quite interesting byproducts. For example, one 
can think of obtaining NSPT results directly in 
momentum space. Also the possibility of exploit- 
ing the so-called Fourier acceleration can be taken 
into account. 

2.4. Dealing with different time scales 

By inspecting Eq. (|h]) one in general does not 
expect the characteristic times to be the same for 
the pure gauge and the fermionic contributions. 
Given the overhead imposed by the inclusion of 
fermions it is important not to waste computing 
time. This is quite easy to do because in order 
to correct for the finite time step e in the integra- 
tion of the Langevin equation we adhere to the 
simplest recipe, that is Euler scheme plus extrap- 
olation e — > 0. Now a first order prescription for 
the integration of an equation of the type 

f(t) = g{t) + h(t) (20) 

1 Notice by the way that the reason why we have dif- 
ferred till now the implementation of fermions in NSPT 
is just the fact that an FFT is easier on APEmille than 
on APE100 



4 



in which different characteristic times are present 
is very easy to work out, a trivial example being 

f(t + e) = f(t) + eg(t) (21) 
f(t + 2e) = f(t + e) + eg(t + e) + 2eh(t + e). 

Having said that, the building blocks represen- 
tation for our algorithm is basically made out of 
three modules: 

• Evolution by the pure gauge contribution to 
the drift F with a certain time step e. 

• Generation of the £V>; this is the only non- 
local piece of the computation, the non- 
locality being anyway traded for multiple 
applications of an FFT, after which every 
operation is trivial. 

• Evolution by the fermionic contribution to 
the drift F with a certain time step e'. This 
does not present any structural difference 
with respect to the first module. 

Of course the second and the third modules are 
tied together and their balance to the first is fixed 
by the considerations sketched a few lines above. 

3. Conclusions, i.e. our shopping list 

As already said, the implementation of the 
program is in progress, so that our conclusions 
are mainly a shopping list of applications we are 
working on. 

• The highest priority is given to the un- 
quenched extension of the a 3 computation 
of the Lattice Heavy Quark residual mass, 
which is a very important building block in 
the determination of the &-quark mass. The 
quenched result (86.2(0.6)(1.0)) has already 
made quite an impact |6| . Being unquenced 
non-perturbative simulations available, the 
unquenced extension of Perturbation The- 
ory is compelling. 

• The new NSPT record is by now the tenth 
order computation of the basic plaquette 
|?j , which strongly confirmed the IR Renor- 
malon dominance picture. Such a picture 



has a clear forecast for the effect of the in- 
clusion of fermions, which is fixed by the 
change in the /3-function coefficients. It is 
worthwhile to verify this. Notice that as a 
byproduct this calls for a high loop deter- 
mination of the critical hopping parameter, 
which could then be compared to the non- 
perturbative computation. 

• Obviously, applications could be devised for 
coefficients needed in various improvement 
programs. We are at the moment consid- 
ering the possibility of extending one loop 
further the computation of the coefficients 
needed in the twisted-mass formulation of 
Lattice QCD [|. 
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